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The multi-scale entanglement renormalisation ansatz (MERA) is argued to provide a natural 
description for topological states of matter. The case of Kitaev's toric code is analyzed in detail and 
shown to possess a remarkably simple MERA description leading to distillation of the topological 
degrees of freedom at the top of the tensor network. Kitaev states on an infinite lattice are also 
shown to be a fixed point of the RG flow associated with entanglement renormalization. All these 
results generalize to arbitrary quantum double models. 
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Renormalization group (RG) transformations aim to 
obtain an effective description of the large distance be- 
havior of extended systems [H. In the case of a sys- 
tem defined on a lattice, this can be achieved by con- 
structing a sequence of increasingly coarse-grained lat- 
tices {Co, Ci, C2, • • • }, where a single site of lattice Ct 
effectively describes a block of an increasingly large num- 
ber tlt ~ exp(T) of sites in the original lattice 
Real-space RG methods can, in particular, be applied to 
study quantum systems at zero temperature, in which 
case each site of Cr is represented by a Hilbert space /Cr 
3 . There the goal is to identify the local degrees of free- 
dom relevant to the physics of the ground state and to 
retain them in the Hilbert space ICr , whose dimension dr 
must be large enough to describe this physics. A severe 
problem of such approach is that in D > 2 dimensions, 
dr must grow (doubly) exponentially in r Q as a result 
of the accumulation of short-range entanglement at the 
boundary of the block. 

Entanglement renormalization 5] is a novel real-space 
RG transformation recently proposed in order to solve 
the above difficulties. Its defining feature is the use of 
disentanglers prior to the coarse-graining step. These are 
unitary operations, acting on the interface of the blocks 
defined by the RG procedure, that reduce the amount 
of entanglement in the system, see figure [T] A major 
achievement of the approach is that, when applied to a 
large class of ground states in both one 0] and two [6] 
spatial dimensions, the dimension dr is seen not to grow 
with T. A steady dr is made possible by the disentan- 
gling step and has several implications d, 0] . It means 
that, in principle, the resulting RG transformation can 
be iterated indefinitely at a constant computational cost, 
allowing for the exploration of arbitrarily large length 
scales. In addition, the system can be compared with it- 
self at different length scales, and thus we can study RG 
flows in the space of ground state or Hamiltonian cou- 
plings. Finally, a constant dr also leads to an efficient 
representation of the system's ground state in terms of a 
tensor network, the multi-scale entanglement renormal- 
ization ansatz (MERA) Q- 
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FIG. 1: RG transformation based on entanglement renormal- 
ization. In order to build an effective site from a block of 
four sites, we first apply disentanglers between sites of the 
block and surrounding sites. In this way part of the short- 
ranged entanglement between the block and its surroundings 
is removed. Then we coarse-grain the four sites into one by 
means of an isometry that selects the subspace IC' C XI®* 
to be kept. We show the case of a tilted square lattice in 
preparation for the toric code where, in addition, each site 
will contain four qubits. 



At zero temperature, strongly correlated quantum sys- 
tems appear organized in a plethora of phases or orders, 
including local symmetry breaking orders and topologi- 
cal orders [1]. Local symmetry breaking phases are de- 
scribed by a symmetry group and a local order param- 
eter, and they are associated with the physical mecha- 
nism of condensation of point-like objects. Transitions 
between two such phases or orders involve a change in 
the symmetry, as described by Landau's theory. A sim- 
ple picture emerges from the perspective of entanglement 
renormalization 31 ■ under successive iterations of the 
RG transformation, ground states with local symmetry 
breaking order progressively lose their entanglement and 
eventually converge to a trivial fixed point, namely an 
unentangled ground state. On the other hand, critical 
ground states describing transitions between these phases 
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are non-trivial — that is, entangled — fixed points of the 
RG transformation. In either case, the MERA provides 
an efficient, accurate representation of the ground state. 

Topological phases are fundamentally different from lo- 
cal symmetry breaking phases Q . They do not stem from 
(the breakdown of) local group symmetries, but their 
topological order is linked to more complex mathemat- 
ical objects, like tensor categories, topological quantum 
field theory, and quantum groups. Physically, topological 
phases exhibit gapped ground levels with robust degen- 
eracy dependent only on the topology of the underly- 
ing space. This, and the fact that excitations above the 
ground level possess anyonic statistics, boosts the inter- 
est of these phases as scenarios for topological quantum 
information storage and processing. Condensation of 
string-like objects (in the so-called string-net models, see 
[9]) has been proposed as a general mechanism controlling 
topological phases. As may be expected, such profound 
differences are also reflected in the way the ground state 
is entangled. Specifically, the notion of topological en- 
tanglement entropy [To'l (the subleading term in a large- 
perimeter expansion of the entanglement entropy of a sys- 
tem) has arisen as a quantitative measure of the ground 
state entanglement due to topological effects. Systems 
with topological order thus provide an unexplored sce- 
nario for entanglement renormalization techniques. 

The purpose of this Letter is to establish entanglement 
renormalization and the MERA as valid tools also for 
the description and investigation of topological phases 
of matter. For simplicity, we analyze in detail Kitaev's 
toric code fu\ , a four- fold degenerate ground state widely 
discussed in the context of quantum computation and 
closely related to Z2 lattice gauge theory and to the 
simplest of Levin- Wen's models for string-net condensa- 
tion Q. We show the following: [i) a MERA with finite, 
constant dr can represent the toric code exactly; (ii) at 
each iteration of the RG transformation, entanglement 
renormalization factors out local degrees of freedom from 
the lattice, while leaving the topological degrees of free- 
dom untouched; (Hi) the MERA representation of the 
four ground states is identical except in its top tensor, 
which stores the topological degrees of freedom; and (iv) 
in an infinite system, the toric code is the fixed point of 
this RG transformation. All these results also hold for 
more complicated models, such as quantum double lattice 
models, that we discuss in the appendix. We conclude 
that the MERA is naturally fitted to represent states 
with topological order, and the entanglement renormal- 
ization offers a new, useful framework for further studies. 

Following we consider a square lattice A on the 
torus, with spin-1/2 (qubit) degrees of freedom attached 
to each link. The Hamiltonian 

-I- □ 

is a sum of constraint operators associated with vertices 



'-f ' and plaquettes '□,' namely 

A+ = 1[X,, Ba = l[Z,. (2) 

ie+ ieO 

Stabilizers act as a simultaneous spin fiip in all 
four qubits adjacent to a given vertex. Stabilizers B\j 
yield the product of group assignments ±1 at the four 
qubits around a plaquette. All stabilizers commute with 
each other and have eigenvalues ±1. Hamiltonian ^ is 
gapped, and states in the ground level (Kitaev states) are 
simultaneous eigenstates of all A^, Bu with eigenvalue 
-1-1. The degeneracy of the ground level (i.e., the number 
of Kitaev states) depends on the topology of the manifold 
underlying the lattice. If this manifold is a topologically 
nontrivial Riemann surface, information is encoded in 
nontrivial cycles, since operators Jliec & '^here Ca.t 
are nontrivial cycles along bonds of the lattice, commute 
with all stabilizers. Besides, such operators along homo- 
logically equivalent nontrivial cycles Ca , Ca have the same 
action on Kitaev states. Hence, for a torus, two logical 
qubits are encoded in the action of these operators. 




FIG. 2: Elementary moves adding plaquettes and vertices to 
a toric code. Arrows stand for CNOT operations. 

Kitaev states are efficiently written in terms of their 
stabilizers. The stabilizer formalism [T3| also provides 
us with a useful language to analyse the action of oper- 
ators on Kitaev states, and has proved instrumental in 
finding an exact MERA. The key observation to this pur- 
pose is that there exist 'elementary moves' [l3|, minimal 
deformations of the lattice and its Kitaev states, that re- 
spect the topological characteristics of the code. These 
moves consists of addition or removal of faces and ver- 
tices together with qubits, and can be written in terms of 
controUed-NOT (CNOT) operators, whose adjoint action 
has a very simple expression in terms of stabilizers: 

I (g) Z ^ Z (g) Z, Z ® I ^ Z ® I, (3) 
I®X^I(gX, X®I X®X. (4) 

Figure [5] depicts the construction of elementary moves. 
The creation of a face is achieved by introducing a new 
spin in a plaquette. Arrows stand for CNOT operators 
from control qubits (all qubits in one of the semiplaque- 
ttes) to the target qubit (the new qubit, introduced in 
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state |0)). The following transformation of stabilizers 
holds (the new site is denoted as n): 

Z1Z2Z3Z4Z5 I — > Z1Z2Z3Z4Z5, (5) 

Zn I > Z\Z-2,Zm (6) 

which ensures plaquette constraints are obeyed. Simi- 
larly, the two relevant vertex constraints are extended to 
the new qubit. The creation of a new vertex is achieved 
instead by introducing a new qubit in state 1+}. This 
qubit now plays the role of control for CNOTs acting on 
the qubits adjacent to one of the split vertices. Stabilizers 
transform as 

XiXiX^Xi^X^ I — > X\XiXj,X^X^^ (7) 

Xn I — X^X\Xm (8) 

which is again compatible with the code constraints. 
Both final sets of operators are the correct stabilizers 
for the code in the modified lattice (remember that 
X^ — Z^ — I.) Similarly, the two relevant paquette 
constraints are extended to the new qubit. 

These operations can be inverted to decouple qubits in 
states |0) and 1+) from the rest system. The disentanglers 
and isometrics, defining both the RG transformation and 
the MERA for the Kitaev states, are made of several of 
these decoupling moves. We regard the original square 
lattice A, on which the toric code is defined, as a (tilted) 
square lattice Cq where each site contains four qubits. 
Then both disentanglers and isometrics act on blocks of 
four sites of Co as in figure [1] — equivalently, on blocks 
of 16 qubits in A. They consist of a series of CNOTs as 
specified in figures [3] and ID 

Upon applying the RG transformation, we obtain a 
coarse-grained lattice Ci which is locally identical to Co 
and where, by construction, the toric code constraints are 
still satisfied. This is quite remarkable. On the one hand, 
it is the first non-trivial example, in the context of entan- 
glement renormalization, where the RG transformation is 
exact , leading to the first non-trivial model that can 
be exactly described with the MERA. On the other hand, 
if we consider an infinite lattice, the above observation 
implies that Kitaev states are an explicit fixed point of 
the RG flow in the space of ground states, as induced by 
the present RG transformation [l^. 

Let us now consider a finite lattice Cq on the torus. 
The coarse-grained state carries exactly the same topo- 
logical information (values of ]^ Z along nontrivial cy- 
cles) as the original state, since the elementary moves pre- 
serve such information at each intermediate step. That 
is, different Kitaev states are not mixed during the RG 
transformation. By iteration, we obtain a sequence of 
increasingly coarse-grained lattices {Cq, £1, £2, • ■ ' 1 ^t} 
for ever smaller toruses. The top lattice Ct will contain 
only a few qubits. Recall that the MERA is made of all 
the disentanglers and isometrics used in the RG transfor- 
mations, together with a top tensor describing the state 
oi Ct It follows that the MERAs for different states 
of the toric code will contain identical disentanglers and 




FIG. 3: (a) The square lattice A for the toric code, with qubits 
(dots) on the links, is reorganized into a tilted square lattice 
£0 where each site is made of four qubits. The lattice constant 
is doubled (dotted lines dissapear) after the RG transforma- 
tion, which produces a new four-qubit site for lattice jC^ from 
every block of sixteen qubits (the twelve light qubits in the 
block are decoupled in known product states), (b) First step 
of the RG transformation: Disentanglers. Arrows stand for 
simultaneous CNOT operators from control to target qubits. 
Disentanglers act on sixteen-qubit domains overlapping with 
four blocks each (thick dashed line, cf. figure [T]) Four qubits 
per block decouple in state |0). 



isometries, and will only differ in their top tensor, where 
all the topological information is stored. 

All the above results automatically extend to the loop 
model considered by Levin and Wen as the simplest of 
their family of string- net models [9j. Indeed, the toric 
code on a square lattice can be locally transformed, us- 
ing the decoupling moves depicted in figure O into a 
toric code on a triangular lattice, which is equivalent to 
the ground state of the loop model defined on the dual 
(hexagonal) lattice. This local transformation shows that 
the topological order of both models are identical, a fact 
already pointed out in 17] and which can also be un- 
derstood in terms of the projected entangled-pair state 
ansatz (PEPS) 18]. 

Finally, our construction generalizes almost straight- 
fowardly to quantum double models (see, e.g., [HI), both 
for Abelian and non-Abelian groups. This is achieved by 
replacing CNOTs with controlled group multiplication 
operators and by paying due attention to the order of 
the operations (see appendix). 

In conclusion, we have shown that several models with 
topological order can be exactly represented with the 
MERA, where topological degrees of freedom are nat- 
urally isolated in its top tensor. We have also seen that 
such models are fixed points of the RG flow induced by 
entanglement renormalization. Our results are an unam- 
biguous sign that entanglement renormalization and the 
MERA, originally developed to efficiently simulate sys- 
tems with local symmetry-breaking phases, provide also 
a most natural framework to study topological phases. 

Acknovirledgements: We thank J. I. Cirac, A. Ki- 
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FIG. 4: (a)-(c) Second step of the RG transformation: Isome- 
tries. (a) Two qubits per block decouple in state |+). (b) Two 
more qubits per block decouple in state |0). (c) One qubit per 
edge, four per block, decouple in state |+). The isometry also 
traces out the twelve decoupled qubits. (d) State of the sys- 
tem after the RG transformation. 
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FIG. 5: Local mapping between the toric code on a square 
lattice (a) and on a triangular lattice (b). The dual model 
in a honeycomb lattice (displayed for reference) is Levin and 
Wen's loop model. 
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APPENDIX A: EXACT MERA FOR QUANTUM 
DOUBLE MODELS 

Here we generalize the above RG transformation and 
MERA to lattice quantum double models (see 
Local degrees of freedom are associated with oriented 
bonds of a lattice A and identified with the group al- 
gebra of a discrete, in general non-Abelian, group G, 
i.e., the Hilbert space spanned by an orthonormal ba- 



sis {\g), g G G}. A change in the orientation of a bond 
corresponds to the map S : \g) i—>- \g~^)- The Hamil- 
tonian is a sum of mutually commuting projectors over 
vertices and plaquettes, 



D(G) 



(Al) 



where vertex projector A„ acts on edges incoming to ver- 
tex V by simultaneous right multiplication by each group 
element, 



' ' heGt 



(A2) 



right multiplication acts as R{h)\g) — \gh)^ and plaquette 
projector Bp selects configurations where the ordered 
product of group elements taken along an oriented cir- 
cuit Cp around p is the unit element of G, 
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FIG. 6: Elementary moves adding plaquettes and vertices to 
a quantum double model. Local degrees of freedom live in the 
group algebra of a discrete group G. The ancilla is initialised 
in state |e) (unit element of G.) Thick arrows stand for con- 
trolled right-multiplication of the target element by the con- 
trol element. Orientation of the edges plays an important role, 
(a) For plaquette addition, operations must be performed in a 
prescribed order (e.g., after application of the arrows in coun- 
terclockwise order (I, II), the new element becomes g\g2-) (b) 
For vertex addition, the new element is initialised in state |e), 
the equal-weight superposition of all elements in G. Here all 
operations can be performed simultaneously. 

Elementary moves are analogous to their counterparts 
for the toric code. The operations generalising CNOTs 
are controlled multiplications by the control element 
(CMs). Figure [5] shows how to create plaquettes and 
vertices using the controlled right multiplication 



A\h, g) = \h, gh), 



(A4) 



where the first element is the control and the second 
element is the target. To cover the case of different 
bond orientations, we also consider the transformations 
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S = (5 (g) 1)A{S ® 1), C = (1 (g) S)A{1 ® S), and 
1) = (S' S)A{S (E) S); explicitly: 

B\h, g) = \h, gh-'), 
C\h, g) = \h, h-'g), 
D\h,g) = \h,hg). (A5) 




FIG. 7: RG transformation for a quantum double model. Ar- 
rows stand for controlled multiplications (CM) from control 
to target elements. Arrow labels denote the type of CM (see 
equations (|A4|I and (|A5I) ). as well as the order in which they 
are applied within each step. The fiducial lattice orientation 
(horizontal bonds pointing to the right, vertical bonds point- 
ing upwards) is assumed, (a) Disentanglers. Four elements 
per block decouple in state |e). (b)-(c) Isometries. In (b), 
two elements per block decouple in state |e). In (c), another 
two elements per block decouple in state |e), and the lattice 
becomes a doubled square lattice with two elements per edge. 
One element per edge, four per block, then decouple in (d) in 
state |e), completing the MERA ansatz. 



By means of these operations, new edges initialised in 

states |e) and 

V heG 

are incorporated into the code, creating new plaquettes 
and vertices. Of course, the inverse elementary moves 
removing plaquettes and vertices from the code, needed 
for the MERA construction, are in general not identical 
to those adding plaquettes and vertices. Note that op- 
erations leading to plaquette addition (or removal) can- 
not be performed simultaneously for non-Abelian groups, 
since the order of multiplication of the elements is impor- 
tant. 

The RG transformation corresponding to a quantum 
double model associated with group G and defined on a 
square lattice proceeds along the same lines as for the 
toric code, but there are qualitative differences. To fix 
the setting, we work with a fiducial orientation of the 
bonds: horizontal bonds are oriented from left to right 
and vertical bonds are oriented upwards. Then: 

• Operations within a plaquette cannot be performed 
simultaneously and must be applied in a certain or- 
der. Hence, disentanglers must be applied in three 
steps, while isometries demand another step with 
respect to the toric code RG. 

• Which of the controlled operations A, B, C, D is 
needed at each step depends on the bond orienta- 
tions. 

The explicit form of the RG leading to a MERA de- 
scription of the quantum double model is shown in fig- 
ure [71 The basic properties of the toric code MERA 
(bounded causal cone, topological degrees of freedom at 
the top of the tensor network, ER fixed point in the infi- 
nite lattice limit) generalise to the quantum double set- 
ting. 
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